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Anomalous temperature dependence of heat capacity of glassy systems during a cooling-heating 
cycle has remained an ill-understood problem for a long time. Most of the features observed in 
the experimental measurement of the heat capacity of a supercooled liquid are shown here to be 
adequately explained by a general model. The model that we propose is motivated by the success 
of landscape paradigm, and describes /3 relaxation in terms of a collection of two-level systems and 
conceives a relaxation as a f3 relaxation mediated cooperative transition in a double-well. The 
anomalous sharp rise in the heat capacity observed during heating is shown to have a kinetic origin, 
being caused by delayed energy relaxation due to nonequilibrium effects. 



The glass transition region is characterized by both 
thermodynamic and kinetic anomalies [3 Q One ex- 
perimentally finds a sharp rise in the measured heat ca- 
pacity of a liquid during rate heating which follows prior 
cooling at a constant rate, the cycle being well extended 
on either sides of the glass transition region QH,EJ. The 
overshoot of the heat capacity is taken to be the signa- 
ture of a glass to liquid transition. While the details of 
the magnitude of the measured heat capacity vary with 
the cooling rate q c and the heating rate q^, the general 
features remain qualitatively the same over a rather wide 
range of q c and q^. In the glass transition region, one also 
encounters with highly nonexponential relaxation of en- 
thalpy, stress and polarization, which is often described 
by the Kohlrausch- Williams- Watts (KWW) form, and a 
very rapid increase of the shear viscosity of the liquid over 
a narrow temperature range Q- 13- 0] ■ Satisfactory expla- 
nation of these anomalies has eluded us for a long time. 
The glassy dynamics is often considered to be a manifes- 
tation of an underlying phase transition However, 
no consensus has still been reached as to the thermo- 
dynamic versus kinetic origin of the observed anomalies. 
Another well-known method to study dynamics in su- 
percooled liquids and glasses is to measure the frequency 
dependence of heat capacity C v (u>), where the relaxation 
time of the heat capacity is the energy relaxation time. 
As discussed by Birge and Nagel, and by Oxtoby, the fre- 
quency dependent heat capacity could provide the much 
needed connection between the thermodynamical and ki- 
netic anomalies |9l Il0|. Oxtoby has, in fact, derived an 
elegant relationship between the time dependence of heat 
capacity and the time dependence of the Active temper- 
ature [lfj. 

It is a common practice to characterize the nonequi- 
librium state of the liquid encountered in time domain 
experiments by the hctive temperature Tf, which, as de- 
fined by Tool and Eichlin [Tj, is the temperature at 
which the nonequilibrium value of a macroscopic prop- 



erty would equal the equilibrium one. If cooling is con- 
tinued through the supercooled regime, the structural re- 
laxation eventually becomes too slow to be detected on 
the experimental time scale, resulting in a limiting hc- 
tive temperature. The limiting Active temperature Tf 
obtained upon cooling is known to depend on q c , while 
the glass transition temperature T g , as measured exper- 
imentally from a cooling-heating cycle, is dependent on 
both q c and qh, a shift to higher values being observed for 
faster rates. If the rates of cooling and heating are taken 
to be the same, that is q c = qn = q, the dependence of T g 
on (|, as shown elegantly by Moynihan et al. £|, is given 
by 



dlnq 



-Ah*/R, 



(1) 



where Ah* can be interpreted as the activation enthalpy 
for the structural relaxation in effect and R is the uni- 
versal gas constant. As pointed out by Moynihan et at, 
it is important for the validity of the above relationship 
that the material be cooled and reheated not only at the 
same rate but the cycle be extended well beyond the glass 
transition range. Ti is also shown to have an identical 
dependence on q c [5|, which has recently been reproduced 
in some model glassy systems [T^ . 

The theoretical analysis of nonequilibrium heat capac- 
ity is a non-trivial problem and has been addressed in 
great detail by Brawer |13t Il4| , Jackie [l5j , and in recent 
time by Odagaki and coworkers |l6j . The two widely used 
expressions of equilibrium heat capacity at constant vol- 
ume are given by 



C V (T) 



dE(T) 
dT 



and 
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where < (AE(T)) 2 > is the mean square energy fluctu- 
ation at temperature T. As is well known, these two are 
equal in equilibrium ergodic system. However, they need 
not be equal in nonequilibrium system. 
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FIG. 1: A schematic representation of the model under con- 
sideration. The horizontal lines within a well represent differ- 
ent excitation levels. Note that the energy levels are in general 
degenerate, as they correspond to the sum of the energies of 
individual TLSs in the collection. 



In this Letter, we present a theoretical analysis of heat 
capacity of a general model of glassy relaxation and show 
that essentially all the features observed during a cooling- 
heating cycle can be explained satisfactorily. The model 
has been conceived in the spirit of the energy landscape 
concept ^3 113 ■ where one describes the system as an 
ensemble of double- well potentials with a broad distribu- 
tion of barrier heights and asymmetries between the two 
minima for the local structural rearrangements [? ? ? ] . 
The model, based on the framework of (3 organized a pro- 
cess, envisages a (3 process as a transition in a two-level 
system (TLS). The waiting time before a transition can 
occur from the level i, labeled either or 1, is taken to 
be random, and is drawn from the Poissonian probability 
density function given by 



ipi{t) = —exp(-t/n), 

Ti 



0,1, 



(4) 



where Ti is the average time of stay in the level i. If pi 
denotes the canonical equilibrium probability of the level 
i being occupied, detailed balance gives the following re- 
lation 



K{T) 



Pl (T) _ n(T) 
Pop) MT) 



= exp[-e/(k B T)}, (5) 



where K is the equilibrium constant for the two levels 
and 1, which are taken to have energies zero and e, 
respectively, and ks is the Boltzmann constant. In our 
model, an a process is conceived as a cooperative transi- 
tion from one well to another in a double-well, subject to 
the establishment of a certain condition. See Fig. 1 for 
a pictorial representation of the model. Each of the two 
wells, labeled 1 and 2, comprises a collection of N (i = 1 
and 2, respectively ) identical, non-interacting TLSs of 
such kind. For a collection of Ni (i = 1, 2) TLSs, a vari- 
able Q(t), (j — 1,2, ,Ni) is defined, which takes on 



FIG. 2: The heat capacity versus reduced temperature plot 
for the system when subjected to a cooling-heating cycle with 
q = 7.5 x 1CP 5 in reduced units. The inset shows the plot of 
the logarithm of the cooling rate q versus the reciprocal of 
the T g . The slope of the linear fit to the data equals —9.04 in 
appropriate temperature units. 



a value if at the given instant of time t the level of 
the TLS j is occupied and 1 if otherwise. Q(t) is thus an 
occupation variable. The variable Qi(t) (i — 1, 2) is then 
defined as 



Ni 



(6) 



Qi{t), which serves to describe the level of instantaneous 
excitation in a collection of TLSs, is therefore a stochastic 
variable in the discrete integer space [0, Ni] . An a process 
occurs only when all the (3 processes (TLSs) in a well are 
simultaneously excited, i.e., Qi — Ni. There is a finite 
rate of transition k from either wells when this condition 
is satisfied. Within the general framework of the model, 
the double- well becomes asymmetric when N\ ^ A^. 

Note that the double-well model has been widely used 
to describe relaxation in glassy liquids and we do not 
specify the microscopic nature of the two states any fur- 
ther. The main new ingredient we introduce is the con- 
dition that transition between the two wells is controlled 
by the collective variable Qi(t). Note also that the dy- 
namics embodied in Fig. 1 is not the standard chemical 
dynamics in a bistable potential because the microscopic 
energy levels are in general degenerate. One should also 
note that in the energy landscape picture the transition 
from the initial state can occur to a limited number (but 
greater than unity) of final states. 

In Fig. 2, we show the heat capacity versus tempera- 
ture C — T curve obtained for our model. The curve looks 
very close to the ones observed in experiments. Note 
the sharp rise in heat capacity during heating. We now 
briefly describe the calculation procedure. The probabil- 
ity Pi (n; T, t) that the stochastic variable Qi takes on a 
value n at temperature T and time t can be shown to 
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satisfy the following stochastic master equation \21 



dPj(n;T,t) 
dt 



= [(JVi-n + l)/76(T)]i^(n-l;r,i) 

+ [(n+l)/7i(T)]^(n + l;r,t) 

- [{Ni-n)/Tn{T)]Pi{n;T,t) 

- (n/niT^Piin; T,t) - k S nM P t (n; T, t) 
+ k 6 n , Ni±1 5j-,i±i Pj(n;T,t), (7) 

where the '+' and '— ' signs in the indices of the Kronecker 
delta are for i = 1 and 2, respectively. The total energy 
of the system at time t can therefore be given by 



E{T,t) 



./Vi iV 2 

^ P x (n; T, t) {N 2 -N x +n) e+V P 2 (n;T,t) n e, 



n=0 



n=0 



(8) 

where the lowest level of the well 2 is taken to have 
zero energy. The set of equations, given by Eq. (7) 
for n — 0,1, ...,Ni and i = 1,2, is solved numerically 
by the matrix method, where the solution is expanded 
in terms of the eigenvectors and the eigenvalues of the 
transition matrix, and the coefficients of the expansion 
are evaluated from the initial distribution. Once we know 
Pi(n;T,t), we can calculate the heat capacity C, as dis- 
cussed below, from an equation, which is essentially a 
form of Eq. (2) modified to incorporate the nonequilib- 
rium effects. 

The system, when subjected to cooling or heating at 
a constant rate, can be envisaged to undergo a series of 
instantaneous temperature changes, each in discrete step 
of AT in the limit AT — > 0, at time intervals of length 
At, whence q = AT/ At Q. A pictorial representation 
of the temperature control during a cooling process with 
finite AT was given by Jackie yj| • If we consider a time 
interval at the beginning of which the temperature has 
been changed from T to T 1 = T + AT, the waiting time 
t bs before an observation can be made is restricted by 
At. The heat capacity C, measured at a time t b s subse- 
quent to a temperature change from T to T' = T + AT, 
is not stationary in time unless t b s is long enough for the 
equilibrium to be established. The measured heat capac- 
ity ( as is the energy ) then becomes a function of the rate 
of cooling / heating as well, apart from T and t b s ■ The 
dependence of C on q c / qh implies that the measured 
heat capacity of a nonequilibrium state depends on the 
history of the preparation of that state. Here we restrict 
ourselves to the case, where q c = qh = q- We therefore 
calculate C(T,t b s ,q) from the following equation: 



C(T,t obs ,q) = 



E(T + AT,t obs ,q)-E{T,0,g) 
AT 



(9) 

where the energies can be obtained from Eq. (8). In the 
present calculation, we have taken t a b s — At. Through- 
out the cycle the transition rates are assumed to be tuned 
with the heat bath temperature T. 

Fig. 2 is the result of the model calculation where 
N% = 6 and N 2 — 10. Temperature T is through- 
out expressed in reduced units of e/fce with e taken 



to be unity. We have set AT = ±0.0015 in reduced 
units. The correspondence to real units is discussed 
later. The model assumes the presence of an energy bar- 
rier for all the intra-well transitions. If ef be the en- 
ergy barrier to the transition from the level i in a TLS, 
the transition state theory (TST) allows one to write 
t,- (T) = h/{k B T)exp[el/(k B T)}, where h is the Planck 
constant. Here we have set e\ — 8 e. We express time 
also in reduced units, being scaled by ri(Tft), where Th is 
the highest temperature considered in a cooling-heating 
cycle. Note that the cycle starts with the equilibrium 
population distribution at Th- The inter- well transition 
rates are equal and independent of temperature. We have 
taken A: -1 = 0.50 in reduced time units. In order to fur- 
ther explore the merit of the model in reproducing the 
experimental results, we have investigated the cooling / 
heating rate dependence of T g for our model. The latter 
has been taken as the temperature of onset of the heat 
capacity increase as observed during heating [23J. The 
Inq versus 1/T 9 plot, as shown in the inset of Fig. 2, is 
linear with negative slope in complete accordance with 
the experimental observations. The slope gives a mea- 
sure of the energy of activation for the relaxation being 
in operation. 

Fig. 3(a) shows a Tf versus T plot for different cool- 
ing rates, where the fictive temperature is calculated in 
terms of energy. The freezing of structural relaxation 
within the experimental time scale is evident from the 
attainment of a limiting fictive temperature. In the inset 
of Fig. 3(a), we show the plot of the logarithm of the 
cooling rate q versus the reciprocal of the limiting fictive 
temperature Tj 1 obtained on cooling. The linearity of the 
plot with negative slope is again in excellent agreement 
with the experimental results. A plot of dTf jdT versus T 
is displayed in Fig. 3(b) for two different cooling / heat- 
ing rates. This is also in accord with the experimental 
observation. 

So far we have presented the results in reduced units. 
It would be useful to have, at this point, an estimate of 
the real terms for reasonable input. For e/ks — 600 K, 
the cooling and heating rates explored here range from 
0.0085 to 0.35 Ks -1 , while the temperature window we 
have looked into lies between 300 K and 120 K. One 
should note that these rates are of the same order of 
magnitude as practised in experiments. 

In order to probe the origin of the behavior of the calcu- 
lated heat capacity, and also of dTf jdT, during heating, 
we plot the energy E versus the heat bath temperature 
T during a cycle, as shown in the inset of Fig. 3(b). The 
fictive temperature evolves in an identical fashion as the 
energy. Note that the energy or the fictive temperature 
goes down during the initial period of heating before it 
starts increasing. The reason is as follows. The presence 
of an energy barrier for all the intra-well transitions re- 
sults in a slow down of the elementary relaxation rates as 
the system is subjected to rate cooling. The system even- 
tually gets trapped into a nonequilibrium glassy state 
on continued cooling. As one subsequently starts heat- 
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FIG. 3: (a) Plot of the Active temperature Tf versus the 
heat bath temperature T in reduced units for different cooling 
rates: ? = 3.0x 10~ 4 , 7.5 x f(T 5 , 2.0 x 10~ 5 , 7.5 x f(T 6 from 
top to bottom. The dashed line traces the Tf = T line. The 
inset shows the dependence of the limiting Active temperature 
Tf obtained upon cooling on the rate of cooling. The slope of 
the linear fit to the data is —9.13 in appropriate temperature 
units, (b) The dTf/dT versus reduced temperature plot for 
the system with Ni = 6 and A^2 = 10 when subjected to a 
cooling- heating cycles. The solid line is for q = 3 x 10 -4 , and 
the dashed line is for q = 7.5 x 10 -5 , both in reduced units. 
The inset shows the evolution of the energy of the system 
during a cooling-heating cycle with q = 7.5 x 10~ 5 in reduced 
units. 



ing, the rate of the elementary relaxation increases, thus 
causing a delayed energy relaxation. This explanation 
further gains support from the calculated heat capacity 
being negative 

Fig. 4 shows a heat capacity versus temperature plot 
during a cooling-heating cycle for our proposed model ob- 
tained with a different set of parameters so chosen that 
the a relaxation becomes more probable within the ob- 
servation time. A repeat of the non-monotonic pattern 
of the heat capacity, although in much smaller scale, is 
notable at higher temperatures during heating for slow 
cooling and heating. However, there has not been any 
report in the literature, to the best of our knowledge, 
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FIG. 4: The heat capacity versus reduced temperature plot 
for the system with iVi = 3 and N2 = 5 when subjected to 
a cooling-heating cycle with q = 5 x 10~ 6 in reduced units. 
The inset shows the same plot for q — 1.5 x 10 -4 in reduced 
units. The axis labels for the inset, being same as those of 
the main one, are not shown. 



of such an observation made experimentally. For faster 
rates, this repeat pattern vanishes as evident from the 
inset of Fig. 4. It is, therefore, reasonable to attribute 
this to the a relaxation. 

A few comments regarding the present work are in or- 
der: 

(1) The hysteresis in the C versus T plot, and also the 
overshoot of the heat capacity observed during heating, 
become progressively weaker as the cooling and heating 
rates decrease, and eventually vanish for sufficiently slow 
rates. This is again in agreement with the long known 
experimental results. 

(2) While the elementary relaxation rates evolve with 
the heat bath temperature, it is the slow population re- 
laxation that gives rise to the nonequilibrium effects. 

(3) The relaxation of energy to its equilibrium value fol- 
lowing a small temperature jump slows down rapidly as 
the temperature is lowered. This is due to the activated 
dynamics assumed here for transition between the two 
levels in a TLS. From the temperature variation of this 
relaxation time, one can have an estimate of the tempera- 
ture where the system starts falling out of the equilibrium 
while cooling at a given constant rate. 

(4) Our model, while simple and microscopic, is quite 
general and also detailed. Its success in reproducing all 
the aspects of experimental results on heat capacity dur- 
ing the cooling-heating cycle is noteworthy. Another im- 
portant aspect is that we need not invoke any singularity, 
thermodynamic or kinetic, to explain the anomalies. It is 
worth mentioning here that a similar model can describe 
many aspects of nonexponential relaxation observed in 
glassy liquids . We are currently investigating the re- 
lationship of the heat capacity anomaly with the fragility 
of the system [2^ . 
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